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The perturbative pion form factor with Sudakov suppression is re-examined. Tak- 
ing into account the multi-gluon exchange in the law Q 2 regions, we suggest that the 



running coupling constant should be frozen at a s (t = J < >) and y < ky > is the 
average transverse momentum which can be determined by the pionic wave function. 
In addition, we correct the previous calculations about the Sudakov suppression fac- 
tor which plays an important role in the perturbative predictions for the pion form 
factor. 
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It is believed that perturbative QCD (pQCD) can successfully describe exclusive pro- 
cesses at asymptotically large momentum transfers JT]. However, the applicability of the 
pQCD to the pion form factor at the present energy is a matter of controversy ||. Recent 
studies [3-5] on the pion electronmagnetic form factor show that the pQCD contributions 
become self-consistent for momentum transfers at a few GeV range. Li and Sterman Q 
give a modified expression for the pion form factor by taking into account the customarily 
neglected partonic transverse momenta as well as the Sudakov corrections. Jakob and Kroll 
H point out that the dependence of the hadronic wave function on the intrinsic transverse 
momentum should be considered in the perturbative calculation. Sudakov corrections come 
from an infinite summation of higher-order effects associated with the elastic scattering of 
the valence partons. However, because the running coupling constant a s becomes rather 
large with b (the distance between a quark-antiquark pair) increasing in the end-point re- 
gions, a cut-off on a s has to be made to evaluate perturbative contributions and to justify 
the self-consistency of perturbative calculations. In this paper, we re-examine the pertur- 
bative pion form factor with the Sudakov suppression. It is pointed out that a s (t) should 
be frozen as t is smaller than a certain value due to the multi-gluon exchange at low Q 2 . 
We suggest that the frozen point is related to the average transverse momentum which is 
determined by the pionic wave function. In addition, we correct the previous calculations 
about the Sudakov suppression factor which plays an important role in the perturbative 
predictions for the pion form factor. 

Let us begin with a brief review on the derivation of the expression for the pion form 
factor in Ref. M Taking into account the transverse momenta that flow from the wave 
functions through the hard scattering leads to a factorization form with two wave functions 
ipi(xi, kyj corresponding to the external pions, combined with a new hard-scattering func- 
tion Th(xi, X2, Q, k^, kr 2 ), which depends in general on transverse as well as longitudinal 
momenta, 

Fn(Q 2 ) = J dxidx 2 J d 2 k. Tl d 2 k T2 ^i(xi, k Tl , Pi) 
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x Tnix!, x 2 ,Q,k Ti , fj,)ipi(x 2 ,k T2 , P 2 ), (1) 

where Q 2 = 2 Pi ■ P 2 , and [i is the renormalization and factorization scale. 
Through Fourier transformation eq. (1) can be expressed as 

Fk{Q 2 )= I dxidx 2 d J°l d ^ 2 (p(xi, bi, Pi, Li) 
Jo \2t^) \2^) 

x T H (xi,x 2 , Q, hi, b 2 , Li)y{x 2 , b 2 , P 2 , At)- (2) 

In this expression, wave function <p(xi, b i? Pj, jj) takes into account an infinite summation of 
higher-order effects associated with the elastic scattering of the valence partons, which gives 
out the Sudakov suppression to the large-6 and small- x regions. At the same time the in- 
trinsic transverse momentum dependence of the wave function provides another suppression 
to the large-6 regions. 

At the lowest order, T H is given by 

167rCiro; s (/i) 



T H (xi, x 2 , Q, k T J — - , 

x ± x 2 Q 2 + (k Tl + k T J 2 

where Cp is the color factor. The wave function can be expressed as 
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ip(x, b, P, n) = exp 
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where 7 g = —a s /n is the quark anomalous dimension in the axial gauge. s(£,b,Q) is the 
Sudakov exponent factor, which reads 
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where 



^ln^QAv^A)], 6 = ln(6A), 



12 ^ 24 

A« = -, A« = - - -n 2 - -n f + -A ln( V), (6) 
3' 9 3 27 * 3 M V 2 ; ' v ; 

where rif is the number of quark flavors and 7 is the Euler constant. 

It should be noted that there are some mistakes in the coefficients of the fourth and the 
sixth terms in b, Q) given by Refs. [4,6]. We find that the correct coefficients should be 

and in P lace of -4§T and in Refs - Ml- lt is that play 

an important role in the evaluation of the pion form factor. In this paper, we examine the 
effects brought about by these corrections. 

Applying the renormalization group equation to Th and substituting the explicit expres- 
sion for Th , we have the following expression for the pion form factor 

Ftt(Q 2 ) = IQttCf / dx\dx2 I b dba s (t) K (y/xiX^Qb) 
Jo Jo 

x <j>(x 1 ,l/b)<j>(x 2 ,l/b)exp(—S(x 1 ,X2,Q,b,t), (7) 

where 



S(x 1 ,x 2 ,Q,b,t) 



2 t 

(s(xi, b, Q) + s(l - x h b, Q)) - — In 
i=i Pi ~b 
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Radiative corrections in higher orders will bring logarithms of the form ln(£//x) into Th, 
where t is the largest mass scale appearing in Th- Ref. 0] points out that a natural choice 
for /i in Th is [A — t and 

t = max (^xix 2 Q, 1/6) • (9) 

If 6 is small, radiative corrections will be small regardless of the values of x because of 
the small a s . When b is large and XiX 2 Q 2 is small, radiative corrections are still large in 
Th, while if will suppress these regions. But with b increasing, a s becomes rather large 
(for example, a s is large than unity when b is large than 0.5/Aqcd GeV -1 for x\ = 0.01, 
x 2 = 0.01 and Q = 2 GeV, see Fig. 1) and accordingly the perturbative calculation loses its 
self-consistency. Therefore, a cut-off on a s is made to evaluate perturbative contributions 
and to justify the self-consistency of perturbative calculation. That is to say, if 50% of 



the result come from the regions where a s is not very large (say, < 0.7), the perturbative 
calculation can be trusted. 

Strictly speaking, the perturbative predictions to the regions where a s is larger than 
unity are unreliable, although these regions are suppressed. In fact, in the regions of small 
X1X2Q 2 and large b, the multi-gluon exchange is important and the transverse momentum 
intrinsic to the bound state wave-functions flows through all the propagators. To respect 
this point, instead of eq. (9) we suggest that 

t = max (y/xiX2Q, l/bp) , (10) 

where 



1/J< k 2 T > if 1/b < J< k T > 



if 1/b > J< k 2 T > 

' (11) 



where J < k T > is the average transverse momentum. With such a choice, the running 



coupling constant will be frozen at a s (t — J < k T >) when b is large and X1X2Q 2 is small. In 
this way, the perturbative contributions to the pion form factor can be calculated from the 
present energy with a reasonable a s , and it should be emphasized that the average transverse 



momentum y < k T > always associates with the hadronic wave function 

The pion wave function: According to Brodsky-Huang-Lepage prescription |7J , one can 
connect the equal-time wave function in the rest frame and the light-cone wave function by 
equating the off-shell propagator in the two frames. They got the wave function [7,3] at the 
infinite momentum frame from the harmonic oscillator model at the rest frame 01 



ifj (a) (x,k T ) = A exp 



k% + m 2 



(12) 



B/3 2 x(l -x) 

where m=0.289 GeV, (3 = 0.385 GeV, A = 32 GeV -1 are parameters which are adjusted 
by using the constraints derived from 7r — > fiu and tt° —* 77 decay amplitudes: 

dx [ — — %ib(x, ky) = ^r- i (13) 
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/ dxijj(x,k T = 0) = ^-, (14) 

JO Jtt 



where f n = 0.133 GeV is the pion decay constant. 

The mean squared transverse momentum is defined as 



where 



<k 2 T >= j ^dx\k T \ 2 \^(x,k T )\ 2 /P q „ (15) 



P«= I ^dx\^(x,k T )\ 2 (16) 



is the probability of finding qq Fock state in the pion. For ip^(x,kx), < ky >=(0.356 
GeV) 2 . Expressing ip^(x,kx) in the 6-space, we obtain 

4>M(x, l/b) = ^x(l - x) exp (-^rr^)) ex P ("2^(1 - x)b 2 ) . (12') 

In order to fit the experimental data and to suppress the end-point contributions for 
the applicability of pQCD, a model for the pion wave function has been proposed in Refs. 
[10,11] by simply adding a factorizing function (1 — 2x) 2 to ip( a \x, kj>). It leads to a wave 
function |TJ[ 

f 14 + ™ 2 1 (17) 



l(3 2 x(l - x) 



ij {b \x,k T ) = A (1 -2a;) 2 exp 
and 

<f>(V( x , l/b) == j0pz{l - x)(l - 2x) 2 exp (- -^-—^ j exp (-2p 2 x(l - x)b 2 ) . (17') 

In the same way as in eq. (12), the parameters are adjusted to be m=0.342 GeV, /3=0.455 
MeV, A=136 GeV" 1 , and < k 2 T > =(0.343 GeV) 2 . 

Numerical calculations : Numerical evaluations for the pion form factor with 0^ and 0^ 
are plotted in Fig. 2. We can find the perturbative predictions are still smaller than the 
experimental data. It is expected to take into account the contributions from higher orders 
and higher Fock states to reach the data at the intermediate energy. 

To evaluate the effects due to the errors in the s(£,b, Q) expression, we adopt the formal- 
ism of Ref. || in our numerical calculations. That is, we choose t as defined in eq. (9) and 



neglect the evolution of with 1/6. In addition, the same two models of the distribution 
amplitudes in Ref. M are used: (a) the asymptotic wave function |T2] 

r s (x) = ^=x(l-x). (18) 

(b) the Chernyak-Zhitnitsky wave function [|13 



t CZ (x) = ^x(l-x)(l-2xf, (19) 

where N c = 3 is the number of colors and f n = 0.133 GeV the pion decay constant. We find 
that the corrected expression (eq. (5)) increases the evaluation of the pion form factor by a 
factor of about 0.8% for the (p as and about 1.0% for the <p cz at Q = 20A QCD . And the effect 
increases with Q decreasing (reaching about 2.0% for (f) as and 3.0% for (j) cz at Q = IOAqcd)- 
The effects are sizable individually for the fourth and sixth term in the s(£, b, Q) expression, 
but fortunately they cancel each other in the final expression. As a result, the whole effects 
on the pion form factor are small. 

Summary: In this paper, we re-examine the perturbative pion form factor with the Su- 
dakov suppression. It is found that in the previous formalism there are regions where a s is 
larger than unity and the perturbative predictions are still unreliable although these regions 
are suppressed. Thus a cut-off on a s has to be made to guarantee the applicability of the 
pQCD. Observing that in the above regions the multi-gluon exchange is important, we sug- 



gest that the running coupling constant should be frozen at a s (t = J < >) when b is large 



and X1X2Q 2 is small by taking into account the average transverse momentum J < >. In 
this way, the perturbative contributions to the pion form factor can be calculated from the 
present energy with a reasonable a s . The essential point of our scheme is that the running 



coupling constant "frozen" is determined by the average transverse momentum y < k|i > 
which always associates with the hadronic wave function. In addition, we correct the previ- 
ous calculations about the Sudakov suppression factor which plays an important role in the 
perturbative predictions for the pion form factor. 
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FIGURE CAPTIONS 



Fig. 1. The evolution of a s with b for Xl = 0.01, x 2 = 0.01, Q = 2 GeV and A QCD = 100 
MeV. The solid line is evaluated with eq. (9). The dashed line is evaluated with eq. 
(10) for <jM\ 

Fig. 2. The pion form factor with 0^ (solid line) and 0^ (dashed line). 
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